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We investigate a class of relativistic fermion theories in 2 < d < 4 space-time dimensions with 
continuous chiral U(lVf) x U(lVf) symmetry. This includes a number of well-studied models, e.g., of 
Gross-Neveu and Thirring type, in a unified framework. Within the limit of pointlike interactions, 
the RG flow of couplings reveals a network of interacting fixed points, each of which defines a univer¬ 
sality class. A subset of fixed points are “critical fixed points” with one RG relevant direction being 
candidates for critical points of second-order phase transitions. Identifying invariant hyperplanes 
of the RG flow and classifying their attractive/repulsive properties, we find evidence for emergent 
higher chiral symmetries as a function of Nf. For the case of the Thirring model, we discover a new 
critical flavor number that separates the RG stable large-A? regime from an intermediate-iVf regime 
in which symmetry-breaking perturbations become RG relevant. This new critical flavor number 
has to be distinguished from the chiral-critical flavor number, below which the Thirring model is 
expected to allow spontaneous chiral symmetry breaking, and its existence offers a resolution to the 
discrepancy between previous results obtained in the continuum and the lattice Thirring models. 
Moreover, we find indications for a new feature of universality: details of the critical behavior can 
depend on additional “spectator symmetries” that remain intact across the phase transition. Impli¬ 
cations for the physics of interacting fermions on the honeycomb lattice, for which our theory space 
provides a simple model, are given. 


I. INTRODUCTION 

The universal behavior of matter near the points of 
continuous phase transitions constitutes one of the most 
intriguing phenomena in statistical physics. Explain¬ 
ing universality was the great success of renormalization 
group (RG) theory, which has by now become one of our 
basic tools to understand systems with many interacting 
degrees of freedom. Systems near criticality can fall into 
universality classes which are characterized by only a few 
basic properties, independent of the microscopic interac¬ 
tions between the particles. In simple bosonic systems, 
the general characteristics that determine the universal¬ 
ity class are well established: the system’s dimensional¬ 
ity, the symmetry of the order parameter, and the pres¬ 
ence or absence of sufficiently long-ranged interactions. 
The reason for this simplicity for a large class of purely 
bosonic theories lies in the existence of just one critical 
RG fixed point, determining the critical behavior of all 
these theories within the corresponding theory space [1] . 

In systems with fermions, however, the situation can be 
more complex, and the question of the defining properties 
of the universality classes could be more subtle. Some is¬ 
sues that arise in these systems are the following: Is the 
above list exhaustive or can one find two systems with the 
same dimensionality, symmetry of the order parameter, 
and range of interaction, but different critical behavior? 
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Does, for instance, the critical behavior depend on addi¬ 
tional (“spectator”) symmetries that do not take part in 
the symmetry breaking pattern? Can there be more than 
one critical RG fixed point in the same theory space? In 
fermionic systems, different types of interactions often 
lead to various possible pairing mechanisms, which al¬ 
low excitations with new quantum numbers, and—upon 
condensation- lead to new collective modes. These new 
modes can be scalar fields as in the BCS theory of super¬ 
conductivity [2], but other types such as vector modes 
are equally well possible [3]. Furthermore, fermionic 
self-interactions are inherently related to Fierz identities, 
with the help of which we can always rewrite any four- 
fermion term as a linear combination of a different set of 
four-fermion interactions. Focusing on only one particu¬ 
lar interaction channel and neglecting all others within a 
single-channel approximation, as is usually done in mean- 
field approaches, involves an ambiguity which may sub¬ 
stantially affect the validity of the approximation [3, 4]. 
In a RG approach, instead, one should incorporate all 
interaction channels that are invariant under the sym¬ 
metry of a given system, and let the dynamics decide 
which one becomes dominant. At the same time, this 
allows to study a whole class of theories in a particular 
theory space and to investigate their decomposition into 
universality classes. 

In this paper, we present a class of relativistic fermion 
field theories in 2 < d < 4 space-time dimensions, which 
allows to study these and related questions. In par¬ 
ticular, we will investigate the space of theories with 
Lorentz, continuous chiral U(iVf) x TJ(Nf), and a set of 
discrete symmetries, with Nf being the number of four- 
component fermion flavors. This includes the ubiqui- 
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tous 3d Gross-Neveu models [5-13], as well as the 3d 
Thiriing model [3, 14-20], both of which have been used 
as testing grounds to study nonperturbative phenomena 
in strongly-coupled fermion field theories, such as chi¬ 
ral symmetry breaking and nonperturbative renormal- 
izability. Lately, these systems receive revived atten¬ 
tion as effective models describing the physics of con¬ 
densed matter systems that incorporate fermionic exci¬ 
tations with relativistic dispersion relation, in particular 
graphene [21-25] and the surface states of topological in¬ 
sulators [26]. Since interactions in graphene are strong 
[27], the question of possible quantum transitions from 
the semimetallic into different Mott insulating [21, 28] or 
superconducting [29] phases has extensively been inves¬ 
tigated previously. We here take a somewhat different 
viewpoint: Instead of focusing on the various possible in¬ 
frared (IR) phases [30] , we target at the ultraviolet (UV) 
structure of our effective-theory space. With the aid of 
functional RG techniques, we map out the fixed points 
and the accompanying relevant and irrelevant directions. 
A fixed point with exactly one relevant direction corre¬ 
sponds to a possible second-order phase transition whose 
critical behavior it governs. We demonstrate that in the 
present theory space multiple such critical fixed points 
may exist, each one of them defining its own universality 
class. Different theories with identical field content and 
microscopic symmetries may therefore be in the domain 
of attraction of different critical fixed points and thus be 
in different universality classes—even though symmetry¬ 
breaking patterns and critical degrees of freedom may be 
completely the same. 

The theory space we propose also represents a sim¬ 
ple example to study the possibility of dynamical emer¬ 
gence of symmetry due to an IR-attractive RG fixed 
point in a higher-symmetric subspace. Emergent sym¬ 
metry is a well-known phenomenon in various condensed- 
matter systems, which often exhibit rotational symmetry 
at low energy, while the microscopic Hamiltonian explic¬ 
itly breaks (continuous) rotational invariance. For lattice 
regularizations of rotational-invariant field theories, dy¬ 
namic enhancement of the rotation symmetry is essential 
in order to avoid the need for fine tuning. The possibil¬ 
ity that Lorentz symmetry, instead of at very high ener¬ 
gies being explicitly realized, could be emergent as a low- 
energy phenomenon is an old idea [31] and a key require¬ 
ment for Horava’s approach to regularize quantum grav¬ 
ity [32]. It is also believed that the critical points in the 
graphene system are Lorentz symmetric [21, 22]. Close to 
a strong-coupling fixed point in (3+l)-dimensional con¬ 
formal field theory [33] and in (2+l)-dimensional [34] and 
(3+l)-dimensional [35] fermionic lattice models super- 
symmetry could be emergent. Emergent supersymmetry 
would be highly desirable for lattice formulations of su¬ 
persymmetric theories, which inevitably break at least 
part of the supersymmetry on the microscopic level [36] . 
The emergence of enhanced internal symmetries has also 
been observed near a fermionic multicritical point with 
S 2 x 0(2) symmetry [37] and in bosonic O(IVi) + 0(N 2 ) 


models [38], leading to the seemingly paradoxical pos¬ 
sibility that Goldstone modes can arise in models with 
initially only discrete symmetries. 

The U (Nf) x U(IVf) theory space we consider in this 
work includes the higher-symmetric U(2Aq) subspace, 
which is the theory space of the continuum Thirring 
model [3, 15]. Already for this simplistic example, the 
question whether or not lower-symmetric perturbations 
out of this subspace are relevant in the sense of the RG 
appears to be nontrivial. It turns out that it can actually 
depend on the number of flavors A): Within our approx¬ 
imation we find that for Nf > 6 the U(2./Vf)-symmetric 
subspace is IR attractive, in accordance with the large- 
Nf analysis [16, 17], while it becomes IR repulsive for 
2 < Nf < 6. Due to an additional Fierz identity the 
Nf = 1 case is special, and we again find that perturba¬ 
tions out of the U(27Vf) subspace are irrelevant, as in the 
large- Nf case. These findings shed new light on previous 
simulation results that employ lattice formulations that 
generically break parts of the microscopic symmetries of 
the continuum theories [13, 16-19]. 

In order to determine the RG flow of our theory space 
we use Wetterich’s functional RG equation [39]. For a 
first analysis, we confine ourselves to the study of point¬ 
like fermionic interactions which is similar in spirit to the 
quantitatively successful derivative expansion for bosonic 
theories. As we shall show, this approximation is equiva¬ 
lent to the usual Wilsonian one-loop RG and as such, 
will become exact to first order in d = 2 + e for all 
Nf, or in any dimension 2 < d < 4 for large Nf. For 
the physically interesting case of small Nf directly in 
d = 3, our simple approximation —as the great major¬ 
ity of all other analytical approaches in the nonpertur¬ 
bative domain—may be not sufficiently controlled. How¬ 
ever, the use of the functional RG in the present case has 
important advantages: first, the method provides mul¬ 
tiple systematic ways to straightforwardly improve the 
present simple approximation used here, e.g, by incor¬ 
porating momentum-dependent vertices or by partial- or 
dynamical-bosonization techniques—all of which are em¬ 
ployed and advanced in recent studies of similar systems 
[3, 40-44]. Second, our approach allows to derive a gen¬ 
eral formula for the one-loop flow of relativistic fermion 
models, which should be directly applicable to systems 
with an arbitrary number of interaction channels. We 
believe that this will be of relevance to future investiga¬ 
tion of more complex systems with critical fermion in¬ 
teractions, e.g., in order to derive an effective theory for 
electrons in the single- [22] or bilayer [45] graphene in 
2+1 dimensions, or in the quantum critical systems with 
quadratic band touching in 3+1 dimensions [46, 47]. 

The rest of the article is organized as follows: In the 
following section we define the theory space that we con¬ 
sider and discuss its symmetries. We briefly introduce 
our method in Sec. Ill which we will use in Sec. IV to de¬ 
rive a general formula for the one-loop beta function of 
relativistic fermion systems with pointlike four-fermion 
interactions. We discuss the flow equations for our sys- 
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tem in Sec. V. Section VI is devoted to general properties 
of the one-loop flow in four-fermion models. In Secs. VII 
and VIII we discuss the fixed-point structure for Nf = 1 
and Nf > 2, respectively. We give an outlook on possible 
phase transitions and critical behavior in Sec. IX, and 
conclude in Sec. X. 


II. FERMION MODELS WITH U(7V f ) x U(7V f ) 
SYMMETRY 


We are interested in the theory space of the U(Ay) x 
U (Nf )-symmetric Gross-Neveu model in 2 < d < 4 space- 
time dimensions, which may be defined by the micro¬ 
scopic action [12] 


Sgn = 


d d a 




— ( 0 Q 0 a ) 5 
2N{ VV V > 


, ( 1 ) 


with space-time index g, = 0 ,..., d — 1 and “flavor” in¬ 
dex a = 1,..., JVf. Summation over repeated indices is 
implicitly understood. We use a four-dimensional rep¬ 
resentation of the Clifford algebra { 7 M , 7 „} = 2<5 AH ,1 4 . 
The Dirac conjugate is given by 0 = A possi¬ 

ble application of this model is the system of Nf species 
of fermions on the honeycomb lattice that interact via 
nearest-neighbor interactions. This interaction can be 
parametrized by the four-fermi coupling g if the four com¬ 
ponents of the Dirac spinor 0 a are associated with the 
electron annihilation operators on the two sublattices and 
the two Dirac points of the honeycomb system, with the 
“flavor” a representing the spin projection [ 21 , 22 ], Ac¬ 
cordingly, Nf = 2 for spin-| particles. In d > 2, the 
coupling g has positive mass dimension, reflecting the 
model’s perturbative nonrenormalizability. However, it 
is now well accepted that the existence of an UV stable 
fixed point ensures that in 2 < d < 4 the model is renor- 
malizable nonperturbatively [5]— a fact that can be rein¬ 
terpreted as maybe the simplest example for Weinberg’s 
asymptotic safety scenario [12]. In the honeycomb-lattice 
system the Gross-Neveu fixed point governs the phase 
transition into the charge density wave phase that is ex¬ 
pected for large nearest-neighbor interaction [21]. The 
one-dimensional theory space defined by Eq. (1) is closed 
under the RG probably to any order [7, 8 ], i.e., once we 
start with a microscopic action of the form ( 1 ) no new 
interactions will be generated by RG transformations. 
This fact has been used to construct the fixed-point po¬ 
tential of the Gross-Neveu model at arbitrary order in 
the fermionic field [48]. However, as we shall see, the RG 
closedness of this model is a special property of the sim¬ 
ple scalar-type interaction, and does generically not hold 
in systems with more complex interactions, such as in the 
Thirring model with a vector-type interaction [15]. We 
will also see that although the Gross-Neveu action defines 
a RG invariant subspace of theory space, infinitesimal 
perturbations on the microscopic level may, depending 
on the number of flavors Nf, be RG relevant and drive 
the system away from the simple Gross-Neveu theory. 


In order to investigate the decomposition of the mi¬ 
croscopic theories into universality classes— beyond just 
relying on the universality hypothesis—and to study the 
stability of the models with respect to perturbations, it 
is therefore mandatory to incorporate the RG flow of all 
operators that are invariant under the given symmetries. 
This defines our U(Xf) x U(ATf) theory space: the space 
of all fermion theories which enjoy (at least) the symme¬ 
tries of the Gross-Neveu model [Eq. (1)]. These are the 
following: 

Relativistic invariance: 

0(a") 1 y 0(a;) 1-4 0 (, (2) 

where 7 MJ/ = §[ 7^,70 and e^ is an antisymmetric 
tensor defining the rotation axis and angle in (2+1)- 
dimensional space-time. Here and in the following, 
where unambiguous, we suppress the flavor index a. 
x' = A0}:r„ with A ^ being the generator of the space- 
time rotations. While the honeycomb lattice explicitly 
breaks Lorentz symmetry, it is, however, expected to dy¬ 
namically emerge close to one of the critical points in 
graphene [21, 22]. Incidentally, the Lorentz transforma¬ 
tion of Eq. (2) can be understood as a particular sub¬ 
group of the global spin-base transformations SL(4, (D) 
[49] constituted by the similarity transformations of the 
Clifford algebra [50]. 

Flavor symmetry: 

0“ U ab ip b , 0“ 1-4 0 b (t/ t ) ba , (3) 

with the unitary matrix U £ U(lVf). Flavor symmetry is 
generated by the (generalized) Nf x Nf Gell-Mann ma¬ 
trices A i, i = 1,..., Nf — 1, together with the identity 
Ao = IjVf- 

Continuous chiral symmetry: Due to our four¬ 
dimensional reducible representation of the Clifford alge¬ 
bra there now exist two additional Dirac matrices, which 
anticommute with all three 7 + 73 and 75 . Their Her- 
mitean product 735 := *7375 generates the continuous 
chiral symmetry 1 

0 a 1 —y e^ 735 0 a , 0 a H 4 0 a e" 1 ® 735 , (4) 

where 9 = 9(a) may depend on the flavor index a = 
1,..., Nf. Continuous chiral symmetry and flavor sym¬ 
metry together are therefore generated by the 2Nf ma¬ 
trices {A 0 ,..., Ajvs-i} <g> { 14 , 735 }, forming the global 
U(lVf) x U(IVf) symmetry of the Gross-Neveu model 
in the reducible representation of the Clifford algebra. 
There is an alternative way to understand this symme¬ 
try [23]: By making use of the orthogonal projectors 
Pl/r ■= }(1 ± 735 ) with Pl /R = P L/R , PlPr = 0, and 


1 Our Dirac matrix conventions are identical to those of, e.g., [3, 
15] with the identification, i.e., simple renaming, 73 «-»■ 74 and 
735 «“>■ 745- 
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Pl + Pr = 1 , we may decompose the four-component 
spinor if) into left- and right-handed Weyl spinors 'j/'l/r : = 
-Pl/rV' and V’l/r := V’Pl/Rj each now representing two 
fermionic degrees of freedom. [On the honeycomb lattice 
" 0 L (V'r) represents the quasiparticle excitations near the 
left (right) Dirac cone.] The Gross-Neveu action then 
decomposes into the two independent parts 


Pgn 
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2Nf 



+ (L R), 
(5) 


without any mixing between left- and right-handed 
spinors. (Note that, in contrast to the usual definition of 
the Weyl spinors in four dimensions, the chiral projector 
Pl/r commutes with 70 , such that t/'lV’r = 0.) Eq. (5) 
is evidently invariant under flavor rotations of the Weyl 
spinors: 

V£^V£K) ba , (6) 

V’r l- t PrVri V’r l- t V’r(Pr) 60 ) (7) 

where the two unitary matrices Pl/r € U(JVf) rnay be 
chosen independently. This constitutes the Weyl repre¬ 
sentation of the U(iVf) x U(7Vf) symmetry. By separating 
the trace and the traceless part of the symmetry genera¬ 
tor, we may as well split the symmetry group as 

U(JVf) x U(JVf) ~ SU(JVf) x SU(iVf) x U(l) v x U(1) A , 

( 8 ) 

where the phase rotations U(l)y (axial transformations 
U(1 )a) correspond to the transformation ( 6 )—(7) with 
scalar matrices Pl/r = e z0 l (Pl/r = e ±lS l); the SU(lVf) 
factors are the remaining transformations with traceless 
generators. For the graphene system with Nf = 2, these 
factors have the following physical meaning: U(l)v corre¬ 
sponds to charge conservation, U(1 )a denotes the trans¬ 
lational symmetry on the honeycomb lattice [ 22 ] , and the 
two SU(2) factors correspond to independent spin rota¬ 
tions in the two Dirac-cone sectors [23]. While the two 
U(l) symmetries are expected to hold in an effective low- 
energy theory of the honeycomb lattice system, the latter 
may possibly (for large coupling) not. Microscopically, 
only the single SU(2) transformation rotating the spin 
simultaneously in both sectors is a symmetry. The sec¬ 
ond SU(2) factor is, however, believed to be emergent if 
the on-site interaction does not become too large [21, 23]. 

TL- 2 . chiral symmetry: 

>->• j5i>, !->■ -^ 75 • ( 9 ) 


Note that different definitions of parity symmetry are in 
principle possible, when combining ( 10 ) with chiral trans¬ 
formations (4). On the honeycomb lattice, the above 
form corresponds to the reflection symmetry which ex¬ 
changes the two Dirac points while not exchanging the 
sublattice labels [22], The second reflection symmetry of 
the honeycomb lattice, which exchanges the spinor com¬ 
ponents belonging to the two different sublattices, is ob¬ 
tained by a combination of (9), (10), and the rotational 
symmetry from ( 2 ) with e= 7 r(<J Ml < 5„2 - 8 ^ 28 ^ 1 ). 

Further discrete transformations may be defined, such 
as time reversal and charge conjugation, both of which 
leave the Gross-Neveu action (1) invariant. For simplic¬ 
ity, we do not list them here, since they do not lead to 
any further restraints on possible operators in the the¬ 
ory. The U(iVf) x U(iVf) Gross-Neveu theory space is 
uniquely determined by the above given set of symmetry 
transformations. 

We now classify all pointlike operators up to the four- 
fermion level with respect to their symmetry. Let us start 
with the two-fermion terms, which represent the building 
blocks of the higher-order operators. Flavor symmetry 
ensures the form 


il) a Oip a ( 11 ) 

with a 4 x 4 matrix O. A basis in the 16-dimensional 
space of 4 x 4 operators is given by the gamma matrices 
and their products: 

O G span{l 4 , 7 ^, 73 , 75 , 7 ^, 735 ,^ 73 ,^ 75 } (12) 

The two-fermion term (11) will be invariant under the 
continuous chiral symmetry (4), if [ 0 , 735 ] = 0, which 
restricts O G span{l 4 , 7 ^, 7 ^, 735 }. For invariance 
under discrete chiral symmetry (9), however, we de¬ 
mand O to anticommute with 75 , which requires O G 
spanjq^, 73 , 735 , 17 ^ 75 }. Finally, parity invariance im¬ 
plies the commutation relation [ 0 , 17175 ] = 0 and thus 
O G span{l 4 , 70 , 72 , 73 , 702 , * 7 o 73 , * 7273 , * 7 i 75 }- If we 
additionally require Lorentz invariance, these restrictions 
mutually exclude each other, and there exists no two- 
fermion term that is invariant under all symmetries of the 
U(lVf) x U(./Vf) Gross-Neveu model. If we relax the con¬ 
dition of Lorentz symmetry down to the invariance un¬ 
der just spatial rotations, the only invariant two-fermion 
term would be ipjoip. In the honeycomb-lattice system 
this represents the quasiparticle density, which vanishes 
at half filling. 

On the level of four-fermion terms there are in principle 
two different types possible, 


The analogous discrete chiral symmetry with 73 instead 
of 75 can be obtained by combining (9) with a prior chiral 
transformation (4) with fixed 9 = n/2 [51]. 

Parity symmetry: We define parity transformation by 
inverting one spatial coordinate, x = (*o, aq, £ 2 ) h > x' = 
(x 0 ,-*1,0:2), 

Ip(x) i^i'y 5 ip(x’), ij>{x) !->• ^(x')*7i75- (10) 


{iP a 0^ a ){ijj b Q^ b ) and {jp a 04> b ){ip b QV), (13) 

which we will refer to as having singlet and nonsinglet 
flavor structure, respectively. Not all of these terms, 
however, are independent: With the help of Fierz iden¬ 
tities, we can always rewrite terms of the second type 
(with nonsinglet flavor structure) as a linear combina¬ 
tion of terms of the first type (with singlet flavor struc- 
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ture), see below. To begin with, it thus suffices to de¬ 
termine the invariant terms with singlet flavor structure; 
the nonsinglet-type ones can subsequently be obtained 
by Fierz identities. From the above discussion it is clear, 
that only the terms with 0 = Q will have a chance to 
be invariant under the given symmetries: Each of the 
(anti)commutation conditions, which lead to restrictions 
of the two-fermion terms, [ 0 , 735 ] = 0 (continuous chiral 
symmetry), { 0 , 75 } = 0 (discrete chiral symmetry), and 
[ 0 , 715 ] (parity symmetry), divides the 16-dimensional 
space of 4 x 4 matrices into two equally large 8 dimen¬ 
sional subspaces of matrices which do and do not respec¬ 
tively fulfill the particular (anti) commutation relation. 
Invariance of the four-fermion term (13) requires 0 and 
Q to be in the same subspace for each of the symmetries. 
Together with Lorentz invariance, this is not simultane¬ 
ously achievable for all three above symmetries if 0 ^ Q. 
Furthermore, for the flavor singlet term in Eq. (13) to be 
invariant under the continuous chiral symmetry we need 
0 = Q to commute with 735 . A basis of flavor-singlet 
four-fermion terms invariant under the above set of sym¬ 
metries is therefore 


(S ) 2 

: = {rrf, 

(P ) 2 

:= (<T735*n 2 , 

(14) 

(V ) 2 

:= (^ a 'y fl ip a ) 2 , 

CO 2 

:= \{ri^r) 2 - 

(15) 


Our theory space includes several previously investigated 
systems: 

(1) The Gross-Neveu model in the four-dimensional re¬ 
ducible representation of the Clifford algebra (“re¬ 
ducible Gross-Neveu model”) with Lagrangian £ = 

+ gs(S ) 2 has been discussed in Ref. [ 12 ]. 

(2) If we choose a gamma-matrix basis in which 735 = 
( l2 -i 2 ) *he iVf four-component Dirac spinors ip a , 
a = 1,..., N{, can be reduced to 2iVf two-component 
Dirac spinors %®, i = 1,..., 2N{ by means of 

and $ = (_ x ^+iV f ) • (16) 

X and x are related to the Weyl spinors by V'l iV’S. = 

( ±x a+Ni ) and ^ ± 4 ’r = ( =Fx“ +Jv f ) • In terms of 
the 2N{ flavors of two-component spinors \ l the in¬ 
teraction channel (P ) 2 can rewritten as a (standard) 
Gross-Neveu interaction [3], 

(■Pf = (x*x) 2 - ( 17 ) 

The Lagrangian £ = + gp(P ) 2 thus de¬ 

scribes the Gross-Neveu model in the irreducible 
representation of the Clifford algebra (“irreducible 
Gross-Neveu model”). Its critical behavior has been 
studied in Refs. [5—10] . 

(3) The system with Lagrangian £ = rpi'y^d^ip + gv(V ) 2 
is known as the Thirring model and has also been 
subject of several previous investigations [3, 14-20]. 


An analogous discussion for the flavor nonsinglet terms 
[second term in Eq. (13)] is possible, but can be kept short 
with the help of the Fierz identities [15]: 

(S) 2 = -I [( S D ) 2 + ( P D ) 2 + ( V D ) 2 + (A D ) 2 ] , (18) 

( P ) 2 = ~ [( S D ) 2 - ( P D ) 2 + ( y D ) 2 - (A D ) 2 ] , (19) 

(V ) 2 = [3(S D ) 2 - 3(P D ) 2 - (V D ) 2 + (A D ) 2 ] , (20) 

(T) 2 = ~ [3 (S D f + 3 (P D ) 2 - (V D ) 2 - (A D ) 2 ] , (21) 


where we have abbreviated 

(S D ) 2 := (r ^ b ) 2 + (r^ b ) 2 , (22) 

(P D ) 2 := (r^4 > b ) 2 + (r^ b ) 2 , (23) 

{V D f := {n^ h f + \(r^ b ) 2 , (24) 

{A d ) 2 := (^ Q *7^73 ^ b ) 2 + (ip a np~/5^ b ) 2 ■ (25) 


This yields an (invertible) one-to-one correspondence be¬ 
tween vectors in the space of flavor singlet terms with 
their “dual” counterparts in the space of flavor nonsin¬ 
glet terms. Therefore, any invariant flavor nonsinglet 
term must be a linear combination of the dual basis 
vectors ( S D ) 2 , (P D ) 2 , (F D ) 2 , and (A D ) 2 . A full basis 
of fermionic four-point functions in the limit of point¬ 
like interactions is thus given by, e.g., the four terms in 
Eqs. (14)-(15) or the four terms in Eqs. (22)-(25), or a 
combination thereof. Put differently, any four-fermion 
theory with the symmetries of the U(lVf) x U(iVf) Gross- 
Neveu model represents one point in the theory space 
spanned by a set of four basis vectors from Eqs. (14)- 
(15), (22)-(25). In the following, we will investigate the 
RG evolution in this theory space with a particular focus 
on possible fixed points. 

Before proceeding, however, let us make a comment on 
the special case of single fermion flavor (Nf = 1), which 
in the condensed-matter applications corresponds to the 
common simplified model of spinless electrons. In this 
case, we have two further relations between the interac¬ 
tion terms 

(^) 2 = WO 2 + (VW ) 2 = (S ) 2 + (P) 2 , (26) 

{V D ) 2 = (V ^) 2 + \{h»A ) 2 = ( v) 2 + (T) 2 , (27) 

in addition to the Fierz identities. Evaluation of the cor¬ 
responding matrix rank yet shows that only one of these 
two is independent and, e.g., Eq. (26) can be obtained by 
linear combination of Eq. (27) with the Fierz identities 
(18)-(21). Thus, for A) = 1 a “Fierz-complete” basis is 
given by just three independent interaction terms, e.g., 
the flavor singlets ( S ) 2 , (P) 2 , and {V ) 2 . The correspond¬ 
ing relation among the singlet invariants is 

(T ) 2 = -3{S ) 2 - 3(P ) 2 - (V) 2 , for JV f = 1. (28) 

This is in agreement with the previous study of spinless 
fermions on the honeycomb lattice [ 22 ], 
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III. FUNCTIONAL RENORMALIZATION 
GROUP 


In order to compute the RG beta functions, we use the 
functional renormalization group in terms of Wetterich’s 
evolution equation for the effective average action Tk [39] , 


dk^k 



d k Rk(r^ + Rk)- 1 


(29) 


where Tk denotes a scale-dependent effective action as a 
function of an infrared RG scale scale k £ [0, A] (A be¬ 
ing the UV cutoff). This action interpolates between the 
microscopic action S for k —> A, and the full quantum 
effective action T (generating functional of one-particle 
irreducible Green’s functions) for k —>■ 0. The quantity 
q. denotes the corresponding Hessian of the effective 
action, and the function regulator Rk defines the details 
of the regularization procedure. Within the approxima¬ 
tion used in this work, all our results given below turn 
out to be independent of the regularization scheme. For 
reviews on the functional RG, see Refs. [40, 52, 53]. 

While the Wetterich equation (29) is an exact iden¬ 
tity for the effective average action, it is difficult to solve 
it without the use of suitable approximation schemes. 
The four-fermion theories considered in this work ex¬ 
hibit a lower critical dimension of d = 2, at which the 
four-fermion couplings become marginal. In dimension 
d = 2 + e for small e any interacting fixed point will 
therefore be of the order g* = O(e) and thus accessible 
via a (renormalized) perturbative expansion. An analo¬ 
gous argument can be made for a large number of fermion 
flavors Nf. Any higher-order fermionic term, as well as 
momentum-dependent fermionic vertices will then be ir¬ 
relevant at such an interacting fixed point, and the effec¬ 
tive average action in the vicinity of the fixed point will 
have the truncated form 


T k = 


d d 2 


w,r + |f(^ + |^(p) 2 


+ m {v) 2 + m iT)2 


(30) 


where we have used the flavor-singlet basis from 
Eqs. (14)-(15) with scale-dependent four-fermion cou¬ 
plings g s ,k, gp,k , gv,k , and g T .k • For N f = 1, the last term 
oc (T) 2 in Eq. (30) can be rewritten as a linear combina¬ 
tion of the three former terms (S') 2 , (P) 2 , and (V) 2 , and 
Tk would be spanned by just three couplings. We have 
neglected the possibility of a wave-function renormaliza¬ 
tion for the fermions, as the fermion anomalous dimen¬ 
sion is known to vanish in the limit of pointlike fermionic 
interactions [15, 48]. In the following, we will use Eq. (30) 
as an ansatz for T^ to solve the Wetterich equation (29) 
approximately. As all omitted terms are perturbatively 
RG irrelevant, this truncation will become exact to first 
order in e close to the lower critical dimension, as well 
as for large number of flavors Nf. For the physical cases 
with e = 1 and small Nf, Eq. (30) can of course only be 


viewed as the simplest possible truncation within, e.g., a 
systematic expansion of F^ in terms of derivatives. Be¬ 
yond the perturbative domain, higher-order terms (e.g., 
momentum-dependent terms) would in principle have to 
be taken into account, and the stability of our results 
against the inclusion of such terms would have to be ver¬ 
ified. In the framework of the functional RG this can 
be done, e.g., by means of partial [23, 42] or dynami¬ 
cal [3, 43, 44, 54] bosonization techniques, or by various 
decomposition schemes in momentum space [40, 41]. A 
third option is offered by working with full potentials for 
fermion bilinears and determining the solution of the flow 
on all scales on a larger function space including weak so¬ 
lutions [55], as has been used for 3d models in [48]. This 
will be left for future work, and we will here confine our¬ 
selves to the study of the RG flow in the limit of pointlike 
(momentum-independent) four-fermion couplings. 


IV. GENERAL FORMULA FOR 4-FERMI BETA 
FUNCTIONS 


By plugging the ansatz for Tk [Eq. (30)] into the Wet¬ 
terich equation (29), the flow equations for the four- 
fermion couplings can straightforwardly (though possibly 
somewhat tediously) be obtained by equating the coef¬ 
ficients on the left and right hand side of the equation. 
However, it appears worthwhile to cast Eq. (30) into a 
more general form 


r fc 





E 


di.k 

2Nf 


(fOif ) 2 


(31) 


describing a general massless fermion system interact¬ 
ing via several short-range interactions. Let us assume 
that the interactions in Eq. (31) represent a “full basis” 
of fermionic four-point functions in the pointlike limit 
in the sense of Sec. II, thus including all flavor-singlet 
terms which are invariant under a given symmetry. In¬ 
variant flavor-nonsinglet terms can be rewritten as a lin¬ 
ear combination of the above by means of Fierz identi¬ 
ties, and in that sense the basis is “Fierz complete”. As 
we shall see in this section, the functional RG equation 
will allow us to compute the beta functions for this gen¬ 
eral ansatz of four-fermion interactions. The resulting 
formula is readily automatable within a computer alge¬ 
bra system. This will allow to compute the one-loop 
beta functions in various fermion systems in arbitrary 
dimension. E.g., it should be applicable to derive effec¬ 
tive theories for several condensed-matter systems fea¬ 
turing chiral fermions, such as single- [22] or bilayer [45] 
graphene, or in the (3+l)-dimensional systems with lin¬ 
ear [26] or quadratic [46] dispersion relation—systems in 
which generically a large number of interactions are com¬ 
patible with the given microscopic symmetries. 2 


2 We note that while we here have assumed relativistic fermions, 
the crucial point for the derivation of the beta functions is the 









7 


In order to formalize the procedure of equating coeffi¬ 
cients, we introduce the “projector” 


j a,b a, / 3 , 7 ,S 


Xl f ,11 

Slpa H’@ H’s 


tjj=xl )=0 


(32) 


with to-be-defined coefficients Cij. In the functional 
derivatives, the Greek letters a, /3, 7 , 5 = 1,..., c ? 7 refer to 
the spinor index, while Latin letters a, b = 1,... ,N{ de¬ 
note the flavor index. d 1 is the dimension of the represen¬ 
tation of the Clifford algebra, { 7 ^, 7 ,,} = 2<5 M „l t i T • When 
applied to our ansatz for the effective action [Eq. (31)], 
Pi gives 


Ar = 1] Cij Mjig 1 (33) 

3:1 


Here we have split the fluctuation matrix into its field- 
independent propagator part F k = r]j.'u_j ,_ 0 (which 
is straightforwardly invertible) and the fluctuation part 
A+ 2) = T<: 2) — ++ For the four-fermion theories, AT^ 
is quadratic in ^>, ijj, and thus the series in Eq. (38) ter¬ 
minates after n = 2 when applying the “projector” Pi. 
In massless systems also the leading term for n = 1 van¬ 
ishes for reasons of symmetry. The only nonvanishing 
term when plugging Eqs. (37) and (38) into (36) is then 
the quadratic term, 


dkgi = —-Pidk STr 


Ad 2 ) ( r M + i ? fe ) 


-1 


1 2 


(39) 


ensuring that the flow of the four-fermi coupling has the 
usual (one-loop) quadratic form 

^ k 9i = y.-i*’ 1 ) (40) 

j,i 


with the symmetric matrix 

Mjx = N { [FriOjOi)] 2 - TriOjOiOjOi) (no sum). 

(34) 

(Wherever possible, we suppress the scale index k of the 
coupling gi = g^ k in the following to avoid confusion.) If 
we choose the coefficients such that 


^^CijMji = 6 a => = (M 1 ) i j , (35) 

j 


we find that by applying Pi on the Wetterich equation we 
“extract” the beta function for the ith coupling constant, 


dkgi 


Pi(d k T k ) 


1 A STr 
2 


d k R k (T^ + Rk)- 1 


(36) 


This can be further simplified by introducing the scale- 
derivative d k ’■= (d k R k )j^, which acts only on the reg¬ 
ulator’s /e-dependence [52], 

d k Rk{ ri 2) + Rk)- 1 = d k ln(T< 2) + R k ), (37) 


and expanding the logarithm in powers of interactions 


ln(T^ +Rk) = Hr ^ 0 + R k ) 


+ E 


n 


AT ( k \T kfi + r^- 1 


(38) 


fermion propagator’s Dirac-matrix structure ~ 7 ^, which has 
an analogous form in nonrelativistic chiral fermion systems with 
quadratic dispersion relation both in 2+1 [56] as well as 3+1 
dimensions [46]. We believe that an analogous formula as derived 
in this section should also be possible to derive for these systems. 


with the coefficient matrix Aiji to be determined. 
Here, we have used the standard abbreviations Vd = 
\vo[(S d ~ 1 )/(2n) d = ++ 1 7 r d / 2 r(d/ 2 )]- 1 , which arises 
from the angular part of the loop integral, and the dimen¬ 
sionless threshold function +A representing the loop 
integral’s radial part. It incorporates the regulator de¬ 
pendence and can be defined by [15] 


if )d = ~d k I' 

Jo 


„d -4 


dq 2 


[1 + r(q 2 /k 2 )] 


2 ’ 


(41) 


In terms of rescaled dimensionless coupling constants 
gi := 4 the beta function reads 


d t gi 


( d — 2 )<?j + 


1 

TVf 


3,1 


t = ln(fc/A). 

(42) 


By employing the ansatz (31) for T*,, Eq. (39) gives after 
some elementary algebra the coefficient matrix 


Aji = ^^2 c im ^Tr (Oi'Y fl OjO m O a^OjOm) 

m 

- Tr{Oa ll OjO m Oj'y ll O l O m ) 

- Tr{O j O m O l 'y li O j 'y li O l O m ) 

- TA{O l O m Oj lll Oa tl OjO rn ) 

+ N f [- Tr {OapOjOm) Tv{O ail OjO m ) 

+ TriOa^OjOn) Tr(CVy M 0,0 m ) 

+ T[(OjO m OiO m ) Tr {^Oj^Oi) 

+ Tr(OjO m ) Tr(O l 'Y fl O j 'y li O l O m ) 

+ Tr (O t O m ) Tr(O j7fl O l7fl OjO m )\ 

- N 2 Tr (OjO m ) TriOiOrn) Tr^O^O,)}. (43) 

Eqs. (42)-(43) represent the main result of this sec¬ 
tion. They are (up to the rescaling) exactly the one-loop 















beta functions as they would have been obtained within 
the standard Wilsonian momentum-shell RG approach. 
In particular due to its simple possible implementation 
within a computer algebra system, we expect this general 
formula to be valuable also beyond the scope of the mod¬ 
els considered in this work. We will use it in the following 
to compute the flow of the four couplings present in the 
U(ATf) x U(JVf) Gross-Neveu theory space [Eq. (30)]. 


V. FLOW EQUATIONS 


A straightforward evaluation of the traces in Eq. (43) 
gives the beta functions for the pointlike four-fermion 
couplings in the four-dimensional Gross-Neveu theory 
space for Nf > 1 


dtgs 


dtgp 


d t gv 


dtgp 


(d - 2 )gs + — [(—4ATf + 2 )g$+ 

+ gs(2gp + 6 gv + 6 gp) + 8 gvgr\, 

(d — 2 )gp + — [(—ANf + 2)gp+ 

+ gp(2gs + 6gv + 6 g T ) + ^gv + 4 St] ■ 


(d~ 2 )g v + — 


4N { + 2 2 

— gv+ 


+ gv ( - g3s + 2 gp + -gp ) + ^ 9s gr 


4N{ 


J 9 

~9p+ 


(d~2)g T + -^ 

+ 9t (~^9s + 2 gp + -gv ) + 7,9sgv 


(44) 

(45) 


(46) 


(47) 


We observe that these flow equations are symmetric un¬ 
der the exchange of gv gr- They generalize several 
previous RG approaches to relativistic fermion systems: 

Setting g P = g v = g T = 0 defines the U(iVf) x U(iVf) 
Gross-Neveu model with four-component Dirac fermions. 
The remaining flow equation 

d t gs = (d- 2)g s - 4A ^ " gf (48) 

agrees with the previous calculation [12]. We also see 
that this subspace is closed under the RG, and gp, gv, 
and gp will not be generated if they all vanish at the 
initial scale. 

The subspace gs = gr = 0 has also been considered 
earlier, and our flow equations agree with the former 
work [15]. In this subspace the system’s 7L 2 x U(lVf) x 
U(iVf) chiral symmetry is elevated to U(27Vf), generated 
by the 41V f 2 matrices {A 0 ,..., A Jv a_ 1 } ® { 14 , 73 , 75 , 73 s}, 
and RG closedness is protected by symmetry. This is 
the theory space of the Thirring model in 2 < d < 4 di¬ 
mensions [3, 15]. In the following, we will refer to it as 
U(2iVf) or Thirring subspace. 

For Nf = 1, when the theory space is just three- 
dimensional, one of the four couplings in Eqs. (44)-(47) 


can be completely eliminated. Using Eq. (28), we find 
gs(S ) 2 + g P (P ) 2 + g v {V ) 2 + g T (T ) 2 = 

(gs ~ 3 gp)(S) 2 + (gp - 3 g T )(P) 2 + (gv - gr)(V) 2 ■ 

(49) 

If we shift the couplings as 

gs - 3 gp >-»• gs, gp - 3 gr gp, gv - gr gv, 

(50) 

the corresponding shifted beta functions for Nf = 1 then 
indeed become independent of the fourth coupling gp, 

dtgs = (d- 2 )gs + 2gs(—gs + gp ~ gv), (51) 
dtgp = (d - 2)g P + 2g P (-g P + gs + 3 gv) 

+ 4 gv(gv - 2 gs), (52) 

dtgv = (d - 2)g v + 2 g v {gv - fa’s + gp) ■ (53) 

We note that Eqs. (51)-(53) are equivalent to the pre¬ 
vious one-loop flow equations for the spinless-fermion 
system on the honeycomb lattice [22]. 3 The flow equa¬ 
tions thus again reflect the fact that for Nf = 1 differ¬ 
ent values for the coupling gp (in terms of the shifted 
couplings) do not correspond to different physical points 
in theory space, but rather can be mapped onto each 
other by means of the Fierz identities; the linearly de¬ 
pendent set of the four elements (S) 2 , (P) 2 , (U) 2 , and 
(T) 2 , which span the theory space, is no longer a basis 
for Nf = 1, but an overcomplete frame, as used, e.g., in 
the theory of signal processing [57]. Of course, the same 
would happen for Nf > 1 if further linearly dependent 
(e.g., flavor-nonsinglet) terms were to be added to the 
effective action. 

Let us express a warning at this point: When count¬ 
ing the number of (physically distinct) fixed points, and 
the number of their accompanying relevant and irrelevant 
directions (as we shall do in the following) it is thus im¬ 
portant to take all identities between the fermionic terms 
into account and reduce the frame to a (linearly indepen¬ 
dent) basis: Otherwise, e.g., critical fixed points, having 
a single relevant direction in the irreducible basis, may 
develop spurious additional relevant (or irrelevant) di¬ 
rections in the redundant directions, preventing them to 
be identified as critical fixed points in the overcomplete 
frame. On the other hand, if by means of a single- or 
few-channel approximation we a priori neglect particular 
(physical) directions in theory space, which in principle 
would be compatible with the symmetry, a fixed point 
with several RG relevant directions could falsely appear 
as critical fixed point in such a truncated description. 
In this work, we therefore advocate the use of a Fierz- 
complete irreducible basis of theory space in order to find 
the number of fixed points and their relevant and irrele¬ 
vant directions. 


! The equivalence can be readily seen by using the coupling trans¬ 
formation g a = -g v , g C i = gs~9v, and g D2 = g P -2g v , where 
9 a,Cl,D 2 are the three couplings used in Ref. [22], assuming the 
full Lorentz symmetry. 
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VI. GENERAL PROPERTIES OF THE 
ONE-LOOP FLOW 

The topology of the flow is determined by the fixed 
points g* of the RG, where all beta functions vanish, 
dtg\ a - = 0. For the fermionic systems in 2 < d < 4 di¬ 
mensions considered here, they separate the theory space 
into a domain of attraction of the Gaussian fixed point 
g* = 0, and the strong-coupling regime, where the flow 
becomes unstable, signaled by a divergence of the renor¬ 
malized couplings at finite RG scale. In the vicinity of 
a fixed point, the RG flow can be linearized and is then 
governed by the stability matrix B, : j = d(dtgi) / dgj \ g *. 
Let v 1 be the eigenvectors and — 0 7 the eigenvalues of 
the stability matrix, i.e., Bv 1 = — Q I v I with I = 1,..., 4 
(J = 1,2,3) for Nf > 2 (Nf = 1). In the basis {id} we 
can thus integrate the linearized flow 

v 1 (t) = id(0) exp(—0 7 f), (54) 

and v 1 therewith defines a RG relevant (irrelevant) di¬ 
rection, if the exponent 0 7 > 0 (0 7 < 0). We call a fixed 
point “critical”, if exactly one exponent is positive and 
all others are negative. Such critical fixed points may be 
associated to second-order phase transitions. 

The quadratic form of the one-loop beta functions 
[Eq. (42)] causes the flow to exhibit the following sim¬ 
plifying properties: 

(1) Any nontrivial fixed point g* ^ 0 has a positive crit¬ 
ical exponent 0 = d — 2. The corresponding RG 
relevant direction is given by the fixed-point vector 
g* itself. 

(2) The (straight) line connecting the Gaussian with a 
non-Gaussian fixed point is invariant under the RG. 

(3) A plane that contains the Gaussian fixed point g* 0 = 
0 and three non-Gaussian fixed points g\,g^,g* c ^ 0 
is RG-invariant if g\ g c are pairwise linear indepen¬ 
dent. Similar relations hold for higher-dimensional 
subspaces. 

These properties are readily be shown by making use 
of Eq. (42). The stability matrix has the form Bfj = 
(d - 2 )5ij + 2 9k^i,kj an d thus [58] 

E B V9j =(d~ m + 2 E 9lA, kj g* = -{d - 2 )g*, 

j kj 

(55) 

where in the last step we have made use of the fixed- 
point equation d t g \ a * = 0. This proves (1). Property (2) 
is shown similarly: Let A g* with A £ R parametrize the 
line connecting O with the non-Gaussian fixed point at 
g* ^ 0. Then 

d t (\g*) = (d-2)(l-\)(\g*) (56) 

and is thus parallel to g* itself. Ad (3): Without loss of 
generality we can define a basis with the two fixed-point 



FIG. 1. Schematic fixed-point structure for Nf = 1 (not true 
to scale!). Any non-Gaussian fixed point g* y= 0 has a rele¬ 
vant direction along g* itself. There are three critical fixed 
points A, C, and £ with exactly one relevant direction. The 
fixed points A, B, and C are located in the higher-symmetric 
U(2iVf) space gs = 0, which is invariant under the RG (yel¬ 
low/light gray). A and V [together with the fixed point Q at 
(gs, gv, gp) = (oo,0, oo)] are located in the invariant plane 
gv = 0 (cyan/dark gray). 


vectors g* A B being the first two basis elements and all 
other basis vectors pointing out of this plane. The flow 
out of the plane is thus 


dtgt 


9m> 3=0 


= (d— 2 )g.i + QfA. 


/ . 9j-™-i,jl9h 
3,l< 2 


i > 3. 
(57) 


Plugging the three (linear independent) fixed-point con¬ 
ditions dtg* A B c into Eq. (57) gives A^n = A .^12 = 
22 = 0. The plane containing and C is 

thus RG invariant. This reasoning can be generalized 
for n-dimensional subspaces containing n(n + l)/2 non- 
Gaussian (and suitably located) fixed points. 


VII. FIXED-POINT STRUCTURE FOR N { = 1 

For the sake of clarity, we perform the following fixed- 
point analysis in d = 3 space-time dimensions. Note, 
however, that within our approximation we obtain the 
very same results in general 2 < d < 4 upon the appro¬ 
priate rescaling of fixed-point values j 4 (d — 2)g and 
critical exponents 0 1 (d — 2)0. We start with the 
one-flavor case IVf = 1, for which the additional “Fierz” 
identity renders the theory space three-dimensional. Out 
of the 2 3 = 8 possibly degenerate and complex solutions 
of the fixed-point equations [Eqs. (51)—(53)] we find 7 
real and distinct fixed points O, A, B, C, V , £, and 
T. Their locations are depicted schematically in Fig. 1. 
From the equations for general Nf > 1 [Eqs. (44)—(47)] 
we infer that the 8th solution Q of the fixed-point equa¬ 
tions, present for Nf > 1, diverges for Nf — > 1. Fixed 
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TABLE I. Fixed points for Nf = 1, their locations, and num¬ 
ber of relevant directions. GN: Gross-Neveu. 



(ds-9 v>9p) 

#(B > 0) 


0 

(0,0,0) 

0 

Gaussian 

A 

(0,0,1/2) 

1 

irreducible GN 

B 

(0, 3+ ^,^ +1 ) 

2 


C 

(0,-5^,-^i) 

1 

Thirring 

V 

(1/2,0,0) 

2 

reducible GN 

£ 

/ 3(3 —\/5) 3-V5 S-2-/E\ 

\ 4 ’ 4 ’ 2 ) 

1 


T 

/3(3+v / 5) 3+V5 5+2^/5\ 

V 4 ’ 4 ’ 2 / 

2 


Q 

( 00 , 0 , 00 ) 

3 

diverges for Nf — > 1 


points O , A, B , and C lie in the two-dimensional, and 
thus RG-invariant, subspace with gs = 0 (yellow/light 
gray plane in Fig. 1). This plane is in fact the Thirring 
subspace with U(21Vf) symmetry. The fixed-point struc¬ 
ture of the U(2iVf) theory space has been investigated 
earlier [3, 15, 20], showing that A and C are critical 
fixed points with exactly one relevant direction within 
this plane. A is the fixed point which describes the Gross- 
Neveu model in the irreducible two-dimensional represen¬ 
tation of the Clifford algebra with a four-fermi interaction 
that can be rewritten in terms of two-component (Weyl) 
spinors x as (P) 2 = (V’ 735 ' 0) 2 °c (xx) 2 [Eq. (17)]. The 
critical behavior of A is well-known [5-10]. The fixed 
point C describes the chiral phase transition expected in 
the three-dimensional Thirring model and as such has 
been dubbed Thirring fixed point [3, 15, 20]— although 
we emphasize that C does not lie on the Thirring axis 
with pure ( V ) 2 = (ijj 7 M V ’) 2 interaction, but its fixed-point 
action includes both (P ) 2 and (V ) 2 contributions [15]. 
The critical behavior of C has recently also been ap¬ 
proached [3, 20]. Here, we are able to determine whether 
or not small perturbations out of the U(27Vf)-symmetric 
plane are relevant in the sense of the RG. From Eqs. (51)- 
(53) we find the following critical exponents: 

A: 0 = (1,-2,- 2 ), (58) 

C : © = (l, -3 + \/5,-3-\/5 + 5), (59) 

and therefore both U(27Vf) critical fixed points remain 
critical also in the lower-symmetric U(iVf) xU(lVf) theory 
space. In the vicinity of both A at {g*s,gvi9*p) = (0, 0, |) 

and C at (0, — , — v ^~ 1 ) any U(27Vf)-breaking but 

U(lVf) x U(JVf)-preserving perturbation is thus RG irrel¬ 
evant for Nf = 1 , since the only relevant directions of 
each fixed point are the fixed-point vectors g and g£ 
itself. 

A second RG-invariant plane is given by gv = 0 
(cyan/dark gray plane in Fig. 1), and it contains be¬ 
sides A the non-Gaussian fixed point V at ( 3 s, 3 y, 3 p) = 
(|, 0, 0). V describes a theory with pure () 2 four-fermi 
interaction, i.e., the Gross-Neveu model in the reducible, 
four-dimensional, representation of the Clifford-Algebra 
[12]. The 8 th solution Q of the fixed-point equations that 


diverges for Nf —» 1 also lies in this subplane at ( 3 , 0 , 5 ) 
with g oc (Nf - l) -1 . The line ( gs,gv,gp ) = A(l,0,1) 
is therefore also RG invariant, just as the Gross-Neveu 
axes A(l, 0,0) and A(0,0,1). As has been observed previ¬ 
ously [22] , both A and V have only one relevant direction 
in this plane. However, while perturbations out of this 
plane are irrelevant in the vicinity of A, such is not the 
case in the vicinity of fixed point T>: For this fixed point 
the critical exponents are 

V: 0 = (1,2/3,-2), (60) 


and thus the reducible-Gross-Neveu fixed point V has two 
relevant directions in the \J(Nf) x U(lVf) theory space for 
Nf = 1. Since there is no higher symmetry that could for¬ 
bid perturbations oc gv the fixed point T> is not a critical 
fixed point for Nf = 1 and cannot describe a second- 
order phase transition, as long as only one microscopic 
parameter is tuned. 

There is, however, a third critical fixed point, located 
at 


C / * * * \ 

£ : \9si9vi9p) 


( 3(3 — x/5) 3 — y/5 5-2^5\ 

^ 4 ’ 4 ’ 2 ) ’ 

(61) 


which has the same critical exponents as fixed point C in 
our approximation, 

£: 0 = (l,-3 + \/5,-3\/5 + 5) . (62) 

Small perturbations near fixed point T> drive the flow to 
either C or £, depending on the sign of gv- We have sum¬ 
marized our results for the Nf = 1 fixed-point structure 
in Tab. I. 


VIII. FIXED-POINT STRUCTURE FOR Nf > 2 
A. Collision of fixed points 

For Nf > 2 the additional flow equation renders the 
U(iVf) xU(iVf) theory space four-dimensional, thus gener¬ 
ating a more complex fixed-point structure. Generically, 
one expects 2 4 = 16 possibly complex or degenerate so¬ 
lutions of the fixed-point equations [Eqs. (44)-(47)]. We 
find that the number of real (and therefore physical) fixed 
points actually depends on the flavor number Nf. In the 
small-Vf regime 2 < Nf < N^ with 

V f (1) = + i (9872 - 144V3345) 3 

1 

+ | (l234 + 18V3345) 3 
~ 3.76 

there are 12 distinct real fixed points and two pairs of 
complex conjugate solutions. Above Nf > Ay 1 both 
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TABLE II. Selected fixed points for Nf > 2, their locations, 
and number of relevant directions. At Nf \ Ay 1 1 ~ 3.76 
the three fixed points XL, X, and J collide. X and J merge 
and subsequently become complex for Nf < Nf , and ex¬ 
change stability of one direction with the fixed point XL. At 
Nf /*■ N f 2 - 1 = 6, X collides with the Thirring fixed point C, 
and J collides with 1C, again exchanging roles: for Nf > 6, C 
and 1C become critical fixed points. A and V on the two (ir¬ 
reducible and reducible) Gross-Neveu (GN) axes ~ (■0735V’) 2 
and (0'0) 2 , respectively, are critical fixed points for all Nf > 2. 



{ 9 *Si 9 *Pi 9 vi 9 t) 


#(B > 0) 


A 

(0, a, 0, 0) 


1 

irreducible GN 

B 

(0,61,62,0) 


r 3, Nt <6 

12 , Nf > 6 


C 

(0, ci,c 2 ,0) 


[ 2 , Nt <6 
[ 1 , Nf > 6 

Thirring 

V 

(a, 0,0,0) 


1 

reducible GN 

Q 

( 9 , 9 , 0 , 0 ) 


2 


XL 

{hi,hi,h 2 , h 2 ) 

f 1, Nf < Nf 
[ 2 , Nt > AT f (1) 


X 

(*1, *2, * 3 ,*4) ( 

(1, 

[2, 

Af (1) <N { <6 

Nf > 6 


J 

(* 1 , *2, * 4 ,*3) ( 

[ 1, 
[2, 

N f (1) <N f < 6 

Nf > 6 


1C 

( 0 , Cl, 0 , C 2 ) 

1 

[ 2 , Nf < 6 
[1, N { >6 


C 

(/1, h, I2, h) 


4 



pairs simultaneously become real and we find 16 real fixed 
points. A selection of them is listed in Tab. II. For par¬ 
ticular values of Nf these solutions become degenerate, 
i.e., fixed points approach each other in coupling space as 
a function of Nf and eventually collide. In general, when 
two fixed points collide, two principally different situa¬ 
tions are possible: (1) The fixed points can “run through” 
each other as a function of Nf and exchange roles with 
respect to the RG stability of the axis connecting the 
two fixed points. This phenomenon is well known for the 
Wilson-Fisher fixed point in 0 4 theory as a function of 
space-time dimension d, which collides with the Gaussian 
fixed point for d /* 4 and becomes unstable for d > 4 [1] . 
It is also observed for the multicritical fermionic fixed 
point on graphene’s honeycomb lattice as a function of 
flavor number Nf [59]. (2) The other possible situation 
is that the fixed points “merge” and eventually disap¬ 
pear into the complex plane. Such has been observed in 
many-flavor QCD [60], three-dimensional scalar [61] or 
fermionic [62] QED, and recently also in nonrelativistic 
systems with quadratic band touching [46]. 

In our system, we find two values of Nf for which a 
collision of fixed points occurs. The first one is partic¬ 
ularly interesting from a general viewpoint: In the limit 
Nf \ Nf ^ we in fact find two triples of solution, each 
consisting of three non-Gaussian fixed points that ap¬ 


proach each other in coupling space and eventually collide 
simultaneously at two different locations. Let us discuss 
one of these triples in more detail. It consists of the three 
fixed points XL, X, and J, which for general Nf > Ay ' 1 
are located at 

ft: {9*s,9*p,9v>9t) = {hi,hi,h 2 ,h 2 ), (63) 

Z ■ {9si9*p,9vi9t) = (*i>*2,*3,*4), (64) 

J ■ {9*s,9p,9v,9t) = U,k) , (65) 

where h 1>2 = h ly2 {N { ), . 4 = *i. 4 (^f) (M > M) 

are functions of Nf, which we do not display explicitly 
for reasons of readability. The symmetry between the 
locations of X and J is determined by the gv XX 9 t sym¬ 
metry of the one-loop flow equations. For Nf \ Nf 1 '* we 
find ii >2 -X hi and 13.4 -X h 2 , i.e., all three fixed points 
collide. Two of them (X and J) merge and disappear 
into the complex plane for Nf < Nf 1 *, while the third 
{XL) remains as a real fixed point. As can be read off 
from Tab. II, X and J are critical fixed points with ex¬ 
actly one RG relevant direction for ATf above (but close 
to) Ay 1 * and exchange one stable direction with the third 
fixed point XL, which becomes critical only for Nf < N^* . 
Thus, at this triple collision both phenomena, merging 
and disappearing into complex plane as well as exchang¬ 
ing roles with respect to RG stability, are realized. The 
analogous behavior can be simultaneously observed for 
the second triple of fixed points with 2 respectively 3 rel¬ 
evant directions (these fixed points are not among those 
explicitly detailed in Tab. II). To our knowledge, the 
present XJ(Nf) x U(ATf) system represents the first-known 
example displaying such triple collision. 

The second value of Nf at which we find degenerate 
fixed-point solutions is 

n{ 2) = 6 . ( 66 ) 

Here, two simultaneous collisions of the first kind (a la 

( 2 ) 

Wilson-Fisher in d /* 4) occur: For Nf /* Ay - ' we find 
that fixed point X collides with the Thirring fixed point 

C: (9* S ,9*p,9v,9t) = (0,Ci,c 2 ,0), (67) 

with the coordinates cy .2 = ci, 2 {Nf) < 0 , i.e., i\ = 0 , 
i 2 — Ci, iff — c 2 , and *4 = 0 for Nf = Nf . The Thirring 
fixed point, which lies in the higher-symmetric U(2Ay) 
subspace, has two relevant directions for 2 < Nf < 6 and 
exchanges roles with respect to stability of one of its RG 
relevant directions with fixed point X for Nf = Ay ’ = 6 . 
Only for Nf > Nf"' 1 (and Nf = 1, see above) the Thirring 
fixed point is thus a critical fixed point, with X being 

critical below (but close to) Ay ; and developing a second 

( 2 ) 

relevant direction above AT f . The same behavior can be 
observed for the fixed point J that collides and exchanges 
its role with respect to RG stability with the fixed point 

1C: (9s,9*p,9v>9t) = ( 0 , 0 , 0 , 02 ), 


( 68 ) 








12 


( 2 ) 

simultaneously at Nf = lV f . We note that the simul¬ 
taneous collision of two pairs of fixed points, just as the 
fact that both triples collide at the same as dis¬ 

cussed above, is a consequence of the invariance of the 
flow equations [Eqs. (44)-(47)] under the exchange of gy 
and gx , which ensures that if (g* s , g* P , gy, gt^) solves the 
fixed-point equations, so does (g* Sl 5 p, gy), and both 
(not necessarily distinct) fixed points have the same crit¬ 
ical exponents. We discuss the stability of the U(2iVf) 
subspace in the vicinity of the Thirring fixed point in 
detail in subsection VIIIC. 

Just as in the one-flavor case, the two (irreducible and 
re_ducible) Gross-Neveu axes given by pure (^ 735 i />) 2 and 
(■i/i'0 ) 2 interactions, respectively, each define a RG in¬ 
variant one-dimensional subspace, since they contain the 
non-Gaussian fixed points 

A: (. g* s ,g* P ,g* v ,g* T ) = ( 0 ,a, 0 , 0 ) (69) 

and 

Z>: (ffs,ffP,Sv,Sr) = M,0,0), (70) 

respectively, where a = a(N{) > 0. Now, for Nf > 2, 
both A and T> are critical fixed points with exactly one 
relevant direction—in contrast to the one-flavor case, 
where X) exhibited two relevant directions. All other fixed 
points, not listed in Tab. II, have two or more relevant 
directions for all Nf > 2 . 

Let us briefly summarize the results we have obtained 
so far (cf. Tab. II): For 2 < Nf < N^ ~ 3.76 there are 
three critical fixed points A, V, and XL. A pair X and 
J emerges from the complex plane at the location of XL 
when Nf = N^\ For N^ < Nf < N^ there are then 
four critical points A, V, X , and J . At N^ = 6 the fixed 

points X and J collide with C and /C, respectively. For 
( 2 ) 

Nf > N^ we still have four critical fixed points, which 
however are now A, V , C, and /C. 

B. RG invariant subspaces 

The fixed-point structure can further be elucidated by 
considering the flow in subspaces that are closed un¬ 
der the action of the RG. Besides the RG-closed one¬ 
dimensional lines connecting each non-Gaussian fixed 
point with the Gaussian fixed point, we find three RG in¬ 
variant two-dimensional planes, each consisting of three 
non-Gaussian fixed points and the Gaussian fixed point. 
These are 

(1) the Thirring subspace, defined by 

(gs,gp,gv,gr) = (0,gp,gv,0), 

consisting of the fixed point A, i.e., the non-Gaussian 
fixed point of the Gross-Neveu model in the irre¬ 
ducible representation of the Clifford algebra, B at 
g* = (0, b\ 1 62 , 0), and the Thirring fixed point C, 


(2) the Thirring subspace’s equivalent for gy -H- g?, 
given by 

( gs,gp,gv,gr) = (o,gp,o,gT), 

consisting of fixed point A , the Thirring fixed point’s 
equivalent fixed point /C, and a not further specified 
fixed point at g* = (0, &i, 0, b 2 ) (equivalent of B), 

(3) the Gross-Neveu subspace, with 

(gs,gp,gv,gr) = {gs,gp, 0 , 0 ), 

consisting of the fixed points A and T> of the irre¬ 
ducible and reducible Gross-Neveu models, as well as 
the fixed point Q at g* s = gp, 

(4) a not further specified subspace, in which 

(gs,gp,gv,gr) = ( gs,gs,gv,gv ), 

consisting of the fixed points XL , £, and Q. By means 
of the Fierz identities (18)—(21), we can rewrite the 
interactions in this subspace as a linear combination 
of the two Fierz-transformed (“dual”) interactions 
( S D ) 2 and (V D ) 2 : 

gs{S) 2 + gs(P) 2 + gv(V) 2 + g v (T) 2 

= - 7 ){gs + 3gv)(S D ) 2 - -{gs - gv)(V D ) 2 . (71) 

Subspaces (1) and (3) can be associated with the corre¬ 
sponding one-flavor subspaces studied above. The RG 
flow in subspace (2) is completely equivalent to the 
Thirring subspace (1). The invariant ( S D ) 2 -(V D ) 2 sub¬ 
space (4) is new. We have depicted the RG flow within 
the invariant planes (1), (3), and (4) in Fig. 2 for the case 
of Nf = 2. While the flow in directions orthogonal to an 
invariant plane vanishes per deftnitionem, we emphasize 
that the answer to the interesting question whether or not 
small perturbations out of the plane are RG relevant de¬ 
pends on the location considered on the plane, as well as 
on the flavor number Nf. For instance, in the vicinity of 
the fixed point XL (cf. right panel of Fig. 2), which is a crit¬ 
ical fixed point for 2 < Nf < N^\ small perturbations 
orthogonal to the (5’ D ) 2 -(F £> ) 2 -plane are RG irrelevant 
(relevant) for 2 < Nf < Nf 1 ^ (Nf > A^ 1 ^). By contrast, 
near the fixed point £, which for all Nf > 2 has four rel¬ 
evant directions, perturbations are always relevant. For 
all Nf > 2 and in the vicinity of all four fixed points 0 1 
A, V, and Q, the Gross-Neveu plane gy = gp = 0 (mid¬ 
dle panel of Fig. 2) is stable under small perturbations 
out of this plane. (Note that this was not the case near 
the fixed point V for Nf = 1, see above.) In the following 
we discuss in detail the stability of the higher-symmetric 
U(2iVf) subspace against symmetry-breaking perturba¬ 
tions in the vicinity of the Thirring fixed point C (cf. left 
panel of Fig. 2), as this question has been under some 
debate in the past [3, 16, 17, 19]. 
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FIG. 2. Flow in RG invariant planes for Nf = 2. Arrows point towards the IR. Left: Thirring subspace with higher U(21Vf) 
symmetry for gs = gr = 0 [15]. Middle: Gross-Neveu subspace for g v = g T = 0. Right: subspace defined by the Fierz- 
transformed interactions ( S D ) 2 and ( V D ) 2 , given by gs = gp and gv = gr and parametrized by the couplings gs = —1( gs + 
3 gv) (horizontal axis) and gy = —\(gs — gv) (vertical axis). 


C. (In-)stability of Thirring subspace against 
symmetry-breaking perturbations 

The large-A/ analysis [14] shows the nonperturba- 
tive renormalizability of the three-dimensional Thirring 
model with a single interaction parameter gy, which is 
equivalent to saying that at large Nf the Thirring fixed 
point is critical, at least within the U(27Vf) subspace. In 
the context of the Thirring model’s lattice version, at 
large Nf it has been shown that small perturbations that 
break the U(2iVf) symmetry are RG irrelevant and the 
U(27Vf) is IR attractive, at least in the vicinity of the 
Thirring fixed point [16]. Our RG analysis is consistent 
with the large-lVf behavior, but shows that the IR at¬ 
tractiveness of the Thirring fixed point does not reach 
all the way down to Nf = 2. Instead, there exists a “crit- 
ical” number of flavors Nf below which the Thirring 
subspace becomes IR repulsive and the Thirring fixed 
point develops a second relevant direction. Below Nf , 
U(21Vf)-breaking perturbations are RG relevant and drive 
the flow to a different fixed point with a lower symmetry. 

This is visualized in Fig. 3, which shows the RG flow 
in the plane spanned by the two most relevant directions 
at the Thirring fixed point C, for different number of fla¬ 
vors Nf. Depending on where one starts the RG flow 
in the vicinity of the Thirring fixed point one finds a 
finite region of starting values for which the couplings 
flow to zero (Gaussian fixed point). Outside this region 
there is a runaway flow. At the boundaries of this re¬ 
gion the flow runs into a critical fixed point with exactly 
one relevant direction. Only for Nf > Nf" = 6 is this 
the Thirring fixed point itself. Below Nf, these are ei¬ 
ther the fixed points H and T> for 2 < Nf < Nf 1 ' 1 or 
the fixed points I and V for Nf 1%> < Nf < Nf 2 \ Di- 
rectly at Nf the Thirring fixed point merges with I. 


( 2 ) 

We emphasize that iV f should not be confused with the 
chiral-critical flavor number Nf- below which the three- 
dimensional Thirring model is expected to allow chiral 
symmetry breaking [3, 15, 17]. We believe that these 
two phenomena are unrelated and the respective critical 
flavor numbers will most likely not coincide. 


IX. PROSPECTS ON LONG-RANGE PHYSICS 

Within the present fermionic truncation of the effective 
average action it is generically hard to associate a given 
critical fixed point with a specific symmetry-breaking 
pattern and corresponding continuous phase transition. 
Moreover, the general structure of the beta functions as 
discussed in Sec. VI renders the single positive critical 
exponent 0 = d — 2, corresponding to a correlation- 
length exponent of the associated phase transition of 
v = 1/0 = 1 /(rf — 2), independent of Nf. While in¬ 
deed many fermionic universality classes in d = 3, e.g., 
of the Gross-Neveu-type [5-12, 23, 63, 64] seem to point 
to v « 1 (±20%), the insensitivity of the fermionic- 
truncation prediction to the specific transition clearly 
calls for more elaborate techniques to investigate the RG 
flow in the vicinity of a given fixed point. Within the 
functional RG, this is, for instance, possible by suitable 
partial or dynamical bosonization techniques [43, 44]. 
The nature of the interacting phases expected at large 
coupling can also be investigated by computing the flow 
of the order-parameter susceptibilities [65], or the flows 
of full potentials for fermion bilinears [55]. Here, we 
content ourselves with an outlook on possible symme¬ 
try breakings associated with the critical fixed points we 
have found and leave a more detailed analysis for future 
work. 

Some of the fixed points considered in the present 
work have been discussed earlier. Let us start with the 
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FIG. 3. RG flow in the (RG noninvariant!) plane spanned by the Thirring fixed point C and the corresponding eigendirections 
vi ,2 with the largest corresponding exponents 0i,2- This plane is given by g = gj + av i + bv 2 with coordinates a, 6 € 1R and 
basis vectors V 2 oc g J (corresponding exponent 02 = 1) and vi, the latter being the second relevant direction for Nf < 6, 
becoming marginal for N{ = 6 , and eventually irrelevant for Nf > 6. The vertical axis a = 0 thus represents the Thirring 
subspace, while finite a ^ 0 corresponds to U(2jVf)-symmetry-breaking perturbations. By varying the starting points of the 
RG flow in the vicinity of C there is a certain finite region for which the flow runs into the Gaussian fixed point O. Outside this 
region we find a runaway flow. At the boundaries the flow runs into a critical fixed point with exactly one relevant direction. 
This is the fixed point H (T>) if a < 0 (a > 0) for 2 < Nf < Nf ^ (top left and middle), the fixed point X (V) if a < 0 (a > 0) 
for Nf ^ < Nf < Nf 2> = 6 (top right and bottom left), and eventually the Thirring fixed C itself (independent of |a| <C 1) for 
Nf > Nf 2 \ which only then becomes critical (bottom right). Directly at JVf = 6 the Thirring fixed point merges with fixed 
point X (bottom middle). In general, V , H, and X do not lie in the present plane—however, we can observe their “projections” 
onto this plane at the points where the flow in directions parallel to this plane vanishes. These projections are marked by X>', 
H 1 , and X 1 . 


Thirring fixed point C. In the U(2iVf) subspace, the 
Thirring fixed point is always a critical fixed point [15]. 
However, only for Nf = 1 and Nf > Nf ’ 2 ' 1 = 6 + 0(d — 2) 
symmetry-breaking perturbations are RG irrelevant. For 
Nf = 1 the three-dimensional Thirring model is expected 
to exhibit a spontaneous breaking of the “chiral” U(2iVf) 
symmetry at large coupling — gy > —gy c ~ —C 2 (Nf) > 0 
with order parameter (ipip)- The critical behavior has 
been discussed within a functional RG approach [20], 
yielding the following predictions for correlation-length 
critical exponent v and order-parameter anomalous di¬ 


mension 77 

C(Nf = 1) : v k, 1.9 T) k, 1.0. (72) 

For Nf > 6 , where the Thirring fixed point again be¬ 
comes critical also in the presence of U(2IVf)-breaking 
interactions, the ordering presumably no longer breaks 
the chiral symmetry [3], but the exact type of ordering 
is unknown. The ordering is also not clearly identifiable 
for the Thirring fixed point’s equivalent for gy -O- gr for 
Nf > 6 , i.e., the fixed point JC. 

By contrast, the fixed point and the associated criti¬ 
cal behavior of the three-dimensional Gross-Neveu model 
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in the irreducible representation of the Clifford algebra 
(A in our notation) is fairly well known [5-10] . As 
expected from the location of the fixed point on the 
(V’735' ( A) 2 = (XXj axis, a nonvanishing order parame¬ 
ter (' 0735 'f/’) = {xx) 7^ 0 occurs for large coupling gp > 
g* P \A = & and thus breaks the parity symmetry while 
leaving the U(2./Vf) symmetry intact. Here again, i/j cor¬ 
responds to Nf flavors of reducible four-component Dirac 
spinors and \ to the corresponding 2iVf two-component 
(Weyl) spinors [see Eq. (17)]. In the system of interact¬ 
ing fermions on graphene’s honeycomb lattice the phase 
transition into the quantum anomalous Hall state that is 
predicted for large next-to-nearest neighbor interactions 
(at least for Nf = 1 [28]) is expected to be governed by 
fixed point A [22] . The critical behavior (for the example 
of Nf = 2) is determined by the exponents 

A{Nf = 2): v « 0.95... 1.04, r) « 0.70... 0.78, (73) 

where the ranges indicate the different predictions ob¬ 
tained by (4 — e) expansion [6], large- Nf expansion [8], 
Monte-Carlo simulations [9], and functional RG [10]. A 
recent overview of the literature results can be found in 
Ref. [23], 

The fixed point V of the Gross-Neveu model in the re¬ 
ducible representation has been considered in Refs. [12, 
22, 63, 64]. V determines the critical behavior of the 
discrete-chiral-symmetry breaking with order parameter 
(ijjtp) which becomes finite for gs > a. On the honey¬ 
comb lattice, it has been ascribed to the transition into 
the charge density wave phase that is expected for large 
nearest-neighbor interaction [21, 22]. The critical behav¬ 
ior for Nf > 2 coincides with the one in the irreducible 
Gross-Neveu model [Eq. (73)], at least to the order that 
the exponents have been computed so far [10, 12, 63, 64] . 4 
However, we would like to emphasize again that T> devel¬ 
ops two relevant directions in the one-flavor case. There 
is therefore a fundamental difference between the irre¬ 
ducible and reducible Gross-Neveu models for Nf = 1: 
As long as only one microscopic parameter of a model, 
e.g., for spinless fermions on the honeycomb lattice, is 
tuned, only A describes a second order phase transition, 
and it can be driven by increasing gp for small gs and 
gv■ By increasing the microscopic coupling gs and keep¬ 
ing gp and gy small and, say, positive we might find 
a phase transition, which, however, is not governed by 
fixed point T> but instead by the critical fixed point £. 

Within the present approximation it is hard to decide 
which order parameter is induced at fixed point £. A 
simplistic approach often applied [40] is given by keeping 


4 In fact, the RG flows of the reducible Gross-Neveu model studied 
in [12, 63, 64] are in principle identical to those of the irreducible 
Gross-Neveu model [10] within the truncations focusing on the 
dynamics of the bosonic order parameter considered so far. In 
view of our results for the larger theory space, the literature re¬ 
sults for the reducible case for Nf = 1 should rather be considered 
as applying more appropriately to the irreducible case. 


track of the “amount of divergence” of the various con¬ 
densation channels, in order to determine in which chan¬ 
nel the couplings diverges “first”. In our one-loop flow, 
the line connecting £ and the Gaussian fixed point O is a 
RG attractive one-dimensional subspace for Nf = 1 and 
positive couplings. E.g., if we start the flow in the UV 
near the reducible-Gross-Neveu axis at g = (gs, gy, gp) 
with 0 < gv, gp ^ gs and gs above but close to the crit¬ 
ical gg = 1/2, the couplings will always run to the 0 £ 
axis before they eventually diverge at a finite RG scale 
to- From the values of the fixed-point couplings [Eq. (61)] 
we find the following ordering of couplings for t —> to 

— -t 3, — -t 1.38, gs,gp,gv -t oo, (74) 

gv gv 

i.e., the coupling gs diverges “fastest”. From this sim¬ 
plistic analysis one might thus speculate that fixed point 
£ governs a condensation in the ( S ) 2 channel with or¬ 
der parameter (i/jtp). That is, with increase of the mi¬ 
croscopic coupling gs the Gross-Neveu model with re¬ 
ducible, four-component, Dirac spinors should exhibit a 
continuous phase transition beyond which the discrete 
TLi chiral symmetry is spontaneously broken, a predic¬ 
tion that is consistent with the mean-field theory for a 
system with large gs(V’V’) 2 interaction and gs > 0 [12]. 
This, of course, requires confirmation beyond the present 
analysis. In any case, however, we believe there is no rea¬ 
son to expect the critical fixed point £ to exhibit the same 
critical behavior as fixed point A. In our approximation, 
the largest critical exponent is 0 = (d— 2 )+ 0 ((d— 2) 2 ) at 
any critical fixed point. Yet, already the second-largest 
exponents that can be associated to the (universal) cor¬ 
rections to scaling do not coincide: at fixed point A we 
obtain 

A(N f = l): u = 2(d~2) + 0((d-2) 2 ) , (75) 

whereas for fixed point £ we get 

£(iV f = l): u = (3- Vh)(d-2) + 0((d-2) 2 ) . 

(76) 

Beyond our approximation, one would expect already the 
leading exponents, for instance v or r], to receive dif¬ 
ferent corrections at the two inequivalent fixed points. 
If true, the irreducible and reducible Gross-Neveu mod¬ 
els as defined by fixed points A and £ thus represent 
an example of two three-dimensional fermion systems 
which both show spontaneous breaking of S 2 symmetry, 
but differ in their corresponding critical behavior. This 
touches a general issue on universality in fermionic sys¬ 
tems: What are the defining properties that determine 
a specific universality class? Our results suggest that in 
fermionic systems the symmetry of the order parameter 
and the dimension and held content of the given sys¬ 
tem does not yet uniquely define the critical behavior. 
Instead, additional “spectator symmetries” that do not 
take part in the symmetry breaking pattern might also 
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play a decisive role. In our case, this is the U(27Vf) ver¬ 
sus Z 2 x U(lVf) x U(lVf) symmetry that discriminates be¬ 
tween the irreducible and reducible Gross-Neveu models. 
We believe that these general questions on universality 
are an interesting direction for future research, and the 
models presented here constitute a suitable playground 
to study them. 

The fixed point XL that becomes critical for 2 < Nf < 
Nf 11 ss 3.8 has the interesting property that for Nf = 3.5 
it is located exactly on the axis at which all four couplings 
coincide, gs = gp = gv = 9 t < 0. By means of the Fierz 
identities the system on this axis can be written with pure 
( S D ) 2 = ( , ijj a 'i/j b ) 2 + { 4 > a 'y 35 ip b ) 2 interaction [cf. Eq. (71)], 
and positive coupling g$ = —2 gs > 0. One would thus 
expect the ( S D ) 2 condensation channel to become critical 
at this fixed point, and this should remain true also for 
other values of Nf not too far from Nf = 3.5. However, 
the critical behavior of this system is unknown to us. 

We are also not able to make any substantial comment 
on the IR behaviors of the fixed points I and J . Both are 
critical fixed points for n \^ < Nf < Nf 2 \ but cannot be 
associated to one of the possible condensation channels 
in an obvious way. 

X. CONCLUSIONS 

Our study of a general class of relativistic fermion the¬ 
ories in 2 < d < 4 space-time dimensions with continuous 
chiral U(lVf) x U (Nf) symmetry has revealed a network of 
RG fixed points. In a unified framework, our description 
includes a number of well-studied models such as vari¬ 
ous versions of the Gross-Neveu model and the Thirring 
model. If persistent also beyond our simple pointlike ap¬ 
proximation, each fixed point defines its own universality 
class thus facilitating both different continuum limits and 
corresponding microscopic “theories” as well as a diver¬ 
sity of possible long-range phenomena. 

We pay particular attention to those fixed points with 
only one RG relevant direction. These are candidates for 
critical points of a second-order quantum phase transi¬ 
tion. It is one of our main results that the nature of these 
critical fixed points in the present model depends on the 
number of flavors. While the fixed point A corresponding 
to the irreducible Gross-Neveu model is a critical fixed 
point for any number of flavors, the fixed point T> for the 
reducible Gross-Neveu model is only a critical fixed point 
for Nf > 2 whereas it is not for Nf — 1. Most interest¬ 
ingly, the Thirring fixed point C is a critical fixed point 
for Nf = 1 and for Nf > 6. While all these results are in 
agreement with the corresponding large-IVf analyses of 
these theories, they indicate that both the microscopic 
as well as the long-range behavior of these systems can 
change drastically as a function of Nf. 

Within the limit of pointlike interactions, we have 
proven a number of properties of invariant subspaces of 
the RG flow. While subspaces of higher symmetry are al¬ 
ways guaranteed to be invariant subspaces, the structure 


of our RG flows allows to define further criteria that are 
not necessarily related to higher symmetries. Subspaces 
with higher symmetry can lead to the phenomenon of 
emergent symmetry in the long-range physics if this sub¬ 
space is not only invariant but also RG attractive towards 
the IR. Again we found that the emergence of higher sym¬ 
metry can be a flavor-number dependent property. For 
instance, the Thirring subspace of higher U(2 Nf) sym- 

( 2 ) 

metry is attractive only for Nf > iV f and Nf = 1 in the 
vicinity of the Thirring fixed point C. 

As our approach is equivalent to the one-loop (2 + e)- 
expansion, at the very least it establishes the existence of 
Nf ’ and its accompanying qualitatively different behav- 
ior for Nf < Nf One should expect that our estimate 
for its value in d = 3 within the present simple trun¬ 
cation is subject to quantitative improvement beyond 
the pointlike limit. Nevertheless, we find it interesting 
that our one-loop result Nf = 6 + 0 (d — 2) appears 
to be near the number Nj^ tlce = 6.6(1) at which the 
IR observables in the lattice simulations show an abrupt 
change [17]. These simulations employ a lattice formu¬ 
lation of fermions (staggered fermions) that microscopi¬ 
cally breaks parts of the U(2./Vf) symmetry. It is there¬ 
fore of decisive importance for the interpretation of the 
simulation results whether these perturbations are rele¬ 
vant in the RG sense or not. Our analysis suggests that 
the answer to this question might in fact depend on the 
flavor number, with a large-iVf regime in which pertur¬ 
bations are irrelevant and an intermediate-W regime in 
which perturbations become relevant. The boundary is 
given by Nf . It should be interesting for future work 
to establish whether indeed the abrupt change found at 
^lattice j n thg simulations signals- instead of an upper 
bound for chiral symmetry breaking—a change in the 
number of relevant directions at the UV Thirring fixed 
point. Such an interpretation would reconcile the seem¬ 
ing disagreement between N^ tlce and the majority of the 
analytical estimates for A r f x c , which appear to be signifi¬ 
cantly lower than Nj a * tlce [3, 66]. It could potentially also 
resolve the contradiction between the critical behaviors 
found in the continuum Thirring model and its lattice 
version, as discussed in [3]. 

For Nf = 1, our theory space describes the system of 
interacting spinless fermions on the honeycomb lattice, a 
simple model for graphene. This model has thoroughly 
been investigated previously [22, 28, 67, 68]. Within our 
RG approach, we rediscover the simple mean-field phase 
diagram which exhibits besides the semimetallic phase at 
weak coupling two gapped phases, which are approached 
by tuning nearest- and next-nearest-neighbor interac¬ 
tions, respectively, above a strong-coupling threshold. 
The appropriate order parameters are ( ipy.’>) (charge 
density wave phase) and (V> 735 'i/') (quantum anomalous 
Hall phase), respectively, each indicating a spontaneous 
breaking of a symmetry. We can associate these 
second-order quantum phase transitions with their cor¬ 
responding critical fixed points. The transition into the 
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quantum anomalous Hall phase is determined by the 
irreducible-Gross-Neveu fixed point A. located on the 
axis with pure ('i/H'ssV’) 1 2 3 4 5 6 7 8 9 10 11 ) as one would naively also expect 
from mean-field theory. A similar expectation, however, 
could fail in the case of the charge density wave transi¬ 
tion. Our results indicate that this transition is not to be 
associated with the reducible-Gross-Neveu fixed point 2?, 
which is located on the axis with pure ( ipi />) 2 interaction 
and for Nf = 1 has two RG relevant directions. Instead, 
the charge density wave transition may be governed by 
fixed point £, at which all short-range interactions be¬ 
come finite. 

Our results also provide evidence for a new feature of 
universality: whereas universality is conventionally con¬ 
sidered as being governed by the symmetry of the order 
parameter, the dimensionality and the number of long- 
range interactions, our theory space appears to contain 
a possible counter-example: for Nf = 1, there are two 
critical fixed points in the theory space (A and £ in 
our notation) which are likely to be associated with a 
second-order phase transition involving the breaking of a 
7Li symmetry. In the irreducible Gross-Neveu model, the 
phase transition describes a breaking of parity symme¬ 
try, whereas the reducible Gross-Neveu model goes along 
with a breaking of a discrete chiral symmetry. In both 
cases, we have a Z 2 order parameter (real scalar field) 
and the same number of massless fermion modes near the 
phase transition. The main conceptual difference of the 


two models lies in the additional spectator symmetries, 
which remain intact across the phase transition. More 
formally, the phase transition is related to two inequiva¬ 
lent RG fixed points lying in different invariant subspaces 
of the full theory space. Quantitatively, we observe that 
the corresponding critical exponents of the fixed points 
differ. Whereas our simple approximation yields differing 
exponents only for the subleading exponents, there is a 
priori no reason why an improved approximation should 
not lead to differences also for the leading relevant expo¬ 
nent. If this is the case for the phase diagram of spinless 
fermions on the honeycomb lattice, this would implicate 
that the two possible phase transitions lie in different 
universality classes. Even though both transitions would 
go along with the breaking of a discrete Z 2 symmetry, 
they should exhibit a different set of critical exponents. 
The recent advances in overcoming the sign problem in 
lattice Monte Carlo simulations [69] may allow to test 
this prediction in the near future. 
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